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Abstract 

Individual cationic site-energies are explicitly determined from molecular 
dynamics simulations of alkali silicate glasses, and the properties and rele- 
vance of this local energetics to ion transport are studied. The absence of 
relaxations on the timescale of ion transport proves the validity of a static 
description of the energy landscape, as it is generally used in hopping models. 
The Coulomb interaction among the cations turns out essential to obtain an 
average energy landscape in agreement with typical simplified hopping mod- 
els. Strong correlations exist both between neighboring sites and between 
different energetic contributions at one site, and they shape essential char- 
acteristics of the energy landscape. A model energy landscape with a single 
vacancy is used to demonstrate why average site-energies, including the full 
Coulomb interaction, are still insufficient to describe the site population of 
ions, or their dynamics. This model explains how the relationship between 
energetics and ion dynamics is weakened, and thus establishes conclusively 
that a hopping picture with static energies fails to capture all the relevant 
information. It is therefore suggested that alternative simplified models of 
ion conduction are needed. 
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1 Introduction 



The transport of mobile modifier ions in silicate glasses below the glass tran- 
sition takes place in a basically fixed environment provided by the glass 
structure. Besides the disordered structure of the network, the Coulomb in- 
teractions among the mobile ions add significantly to the complexity of the 
underlying energy landscape, because their contribution is both long-ranged 
and non-static. 

A hopping picture of the dynamics has been widely accepted,^ based on 
evidence for the existence of well defined ionic sites, from experiments^'^'^'^ 
and also from simulations^i^"^. Basically hopping models describe the glass 
system by a representation of the effective energy landscape experienced by 
the ions. The conformational disorder of the network is typically represented 
by a distribution of site-energies. They are treated as static, relying on the 
separation of timescales between ion dynamics and structural relaxation.— 
Energy values for individual sites and barriers are normally treated as in- 
dependent of each other in realizations by lattice models; see, e.g.,-. This 
common assumption is however not trivial, and it has been dropped in at 
least one case.~ Numerical investigations of such simplified models have es- 
tablished that the Coulomb interaction between the mobile ions is required 
for a completely correct representation of the dynamics,— and a quantitative 
model of the ionic conductivity has been built with the interactions among 
the ions as the main element.— Most other approaches however disregard 
the Coulomb interaction among the mobile ions.-'^»^>^»^>^»^>^ This imme- 
diately raises the question whether the Coulomb interaction can be simply 
taken into account by a local modification of the energy landscape. 

Recently a statistical method allowed us the identification of all ionic sites 
in realistic molecular dynamics (MD) simulations.— Most of the conclusions 
have been confirmed by other groups.—"^. By direct observation of the ion 
hopping dynamics, key predictions about the ion dynamics could be tested.—. 
Furthermore, different spatial aspects of the conduction paths have been 

Most of the hopping models mentioned above use the concept of site- 
energies. However, typically the properties of these site-energies are pos- 
tulated in an ad-hoc manner. In the present work we express the proper- 
ties of the energy landscape via explicitly determined site-energies. These 
are obtained based on our (MD) simulations of alkali silicate systems, from 
time-averages of the particle energy of the ion residing in a site. Two key 
questions are addressed: What are the properties of the site-energies for a 
real ion conductor? Do these site-energies indeed determine thermodynamic 
or dynamic properties of the corresponding sites? In particular, since the 
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Coulomb interaction between cations can be expected to be strong, it is not 
evident a priori that it can be taken into account as merely a contribution to 
single-particle or, equivalently, single-site energies. This aspect is discussed 
in detail by examining the correlation of site-energies to thermodynamical 
and dynamical observables. 

2 Method 

In our MD simulations the potential energy is given by pair interactions as 
f/sys = ^/'^Y.ij^iUij- The interaction potential Uij{rij) = qiq^ / {A-nsQrij) - 
Aijr~^+Bij exp{—Cijrij) consists of the Coulomb interaction for point charges 

and of a short range part of Born-Mayer type. The qij reflect ap- 
propriately chosen partial charges^. 

For a finite system, single-particle energies Ui can be computed just like 
f/gys as the sum of all pair contributions Uij involving the particle i, Ui = 
Uij- While Ui = 2?7sys, any change Af/gys to the energy of the system 
resulting from a change of particle i is directly given by AUi. Considering a 
cation as particle i, the contribution to Ui only from interactions with e.g. 
other cations can furthermore be isolated simply by restricting the sum to 
the appropriate j. 

The situation is complicated by the use of periodic boundary conditions: 
To avoid finite-size effects, the system is treated as an infinite repetition of 
periodic images of the simulation box.— Due to the long range nature of the 
Coulomb interaction, it cannot be computed by simple summation up to a 
finite cutoff in this case. Instead all periodic images of the simulation box 
must be taken into account. We use the Ewald method— for this purpose, 
where the Coulomb energy is broken into three formal contributions: Only 
the real part U^^ is a pair term like U^^; the Fourier part Ufj^ relates to 
the infinite set of periodic images of particle j, and the self-correction Uf^ is 
globally assigned to particle i only. 

Yet a deeper analysis shows that the Coulomb part U^^ of a single-particle 
energy can still be assembled from these terms, and that it can still be divided 
into contributions from groups of like particles. One can notably still separate 
the interaction of a cation with other cations Uf^^ from the interaction with 
network atoms because Uf^ can be traced exclusively to the cation 

part. A constant correction— for the nonzero charge of the partial systems 
under consideration is neglected, because it does not affect the results of this 
study, where energy differences are relevant. 

As the locations and structural properties of the sites are stable on the 
time scale of alkali ion diffusion, the potential energy of cation i is determined 
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by the site s{i) it occupies. A site-energy E{s) can therefore be identified 
by taking the time average over the particle energy of the ions i{s, t) residing 
in s at any given time t, i.e. E(s) = (f/j(s,j)(t))^. Like the particle energies, 
the total value of the site-energy E^°^{s) is separated into the contributions 
£"^'^*(s) and E^'^^{s). Of course, with this construction E{s) cannot be inter- 
preted as a bare site-energy but refiects the site-energy under the condition 
that it is populated by an ion. 

3 Systems 

We analyzed the site-energies for the alkali disilicate systems (Li20) -2(8102) 
and (K2O) -2(8102), named L8— 2 and K8— 2 respectively. These systems have 
already been described in our previous work^^. They contain 270 cations 
among 1215 atoms, at experimental densities Interactions are given by a 
potential by Habasaki et al.— . The systems were propagated with a timestep 
of 2 fs in the canonical ensemble, using a Nose-Hoover thermostat—. Po- 
sitions and energy contributions were stored every 0.1 ps, with energies 
averaged over five values sampled 20 fs apart. Energy data was produced 
for 10 ns of simulations at 850 K. 

Our analysis identifies for these datasets 291 and 276 sites respectively 
in L8— 2 and K8— 2. Ions are in all cases residing in a site during more than 
98% of the times, according to our algorithm that drops excursions.— 

4 Results 

4.1 Properties of average site— energies 

It turns out that the fundamental characteristics of the energy landscape 
are qualitatively identical in both systems. Note that these energies result 
from long time averages. Possible temporal fiuctuations are discussed further 
below. Broad ranges of site-energies refiect the variation in the disordered 
structure. For the distributions can be fitted by Gaussians, with the pa- 
rameters given in Tab. [H The properties of the distributions are compatible 
with earlier results, determined both for sites^, and without reference to ion 

The infiuence of interactions with both cations and network on the shape 
of the site-energy distributions is demonstrated in Fig. [1] for the system 
L8— 2. For each site the cation contribution E'^^^{s) is plotted against the 
corresponding network energy E^'^^{s). Both sets of energies have standard 
deviations of 5 to 6 eV and cover ranges of more than 20 eV. In the system 
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System /i[eV] OdisAeV] 
LS-2 -8.8 0.33 
KS-2 -5.6 0.24 




KS— 2 standard deviations of 3 eV and ranges of 16 eV are found. In all 
cases the significantly narrower distributions for i?*°*(s) are reached because 
a strong and clear anticorrelation between E'^^'^{s) and i?°''*(s), evident in 
Fig. [H reduces the spread of values for the resulting total site-energies. 

Naturally, the strong variation of can be interpreted as a conse- 
quence of the disordered network structure. The presence of non-bridging 
oxygens (NBO) gives rise to a negative partial charge. Thus, due to the 
strong Coulomb interaction, a minor variation in the distance between the 
site and the nearest NBO may result in a significant variation of the network 
energy. The typical distance between a Li-ion and a NBO is 2A. For example, 
variation of this distance by O.lA gives rise to an energy variation of ±0.4 eV, 
for the potential used in these simulations. 

To test the assumption of many lattice models that the site-energies are 
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Figure 2: Correlation among the energies of neighboring sites in system LS— 2 

uncorrelated, the energies of all sites were compared to those of their next 
neighbors, defined by the first nearest-neighbor shell of the alkali pair distri- 
bution. The results are shown in Fig. [2] for the glass LS— 2, whose behavior is 
again typical for all investigated systems. For each site s the average energy 
of the neighbors of s is plotted against the site-energy of s itself. Data for the 
contribution E"^^^ of the matrix alone is given by the open circles. They show 
a strong spatial correlation of the network energies between adjacent sites. 
The same correlation is found for the cation interactions E'^^^, as it may be 
expected when E^^^ and E^'^^ are anticorrelated themselves. The total ener- 
gies are shown as squares in the inset. For them spatial correlations are much 
weaker (regression coefficient 0.26 ± 0.1), because the remaining imperfec- 
tions in the cancelation of E^'^^ and E^^^ introduce sufficient fluctuations into 
^tot Qj-eate a nearly uncorrelated disordered landscape. 

4.2 Temporal variations of site— energies 

Temporal variations of site-energies can occur as a consequence of different 
physical effects: (1) Gradual structural relaxation of the silicate glass, (2) 
Fluctuations due to motion of adjacent ions. Both aspects are analyzed in 
this part. 
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The long-time variation of E{s), resulting from the structural relaxation, 
was estimated by comparing the mean energy values for each site from the 
first and from the second half of the available 10 ns of data. Judging from the 
width of the Gaussian distribution of these energy differences, the residual 
long-time fluctuations of a site-energy are 0.03 eV and 0.02 eV for the 
Lithium and Potassium glasses. These fluctuations, related to systematic 
drifts are small as compared to the overall width of the Gaussian distribu- 
tions; see Tab. I. 

As mentioned above the energy value for a single site results from an 
average over the whole simulation run. Random fluctuations in this site- 
energy can be observed by defining the energy value instead from very short 
averages. Two time scales have to be compared. First, the time scale Tfiud 
characterizes the time where typical energy fluctuations of a site are explored. 
Second, r^es describes typical residence times of an ion in a site. The distri- 
bution of Tres has a median of a few hundred picoseconds (see below). Only 
if '^fiuct Tres it is justified to characterize the properties of a site, as experi- 
enced by an ion, via the average E{s). In the opposite limit one would have 
to introduce a time-dependent energy E{s). 

In order to check how stable and distinguishable the individual site- 
energies are, we plot in Fig. [3] the average growth of the standard deviation 
aE{t) for the set of energy values that are incorporated into the mean site- 
energy. The variable t denotes the width of the interval, used for calculating 
the average. The development of the total energy E^°^ is shown in the upper 
part, the lower part treats the network contribution E^^^. The cation part 
£jcat omitted, because its behavior is very similar to that of E^^^, due to 
the correlation demonstrated above. 

First, it turns out that for LS— 2 beyond 100 ps a nearly constant value of 
<JE(t) is observed. Thus, indeed, the typical fluctuation time Tjiuct is much 
shorter than Tres- The range of values sampled for the energy of a typical 
site has a considerable width. It reaches a value of a{t = 4.6 ns) = 0.31 eV at 
the end of the plot for E[f in LS-2, and of a{t = 4.6 ns) = 0.35 eV for E^^ 
in KS— 2. The overall stability of the energies on a timescale of more than 5 
ns is supported by the small errors for determined above. 

The site-energies, introduced so far, provide an average view of the energy 
landscape, as it is adopted by most hopping models. In principle one might 
envisage a scenario where dynamical events induce systematic changes in the 
local configuration. Specifically the energy of a site might change as a con- 
sequence of an ionic jump. For example, in the MIGRATION concept™ an 
ion, jumping to a site, is thought to be gradually stabilized by the subsequent 
adaption of the other ions to that jump. In this scenario one might expect 
a gradual decrease of the site-energy via a decrease of E'^^^. Of course, for 
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Figure 3: Growth of the standard deviation cje of energies for a site with 
sampling time. The behavior of E'^^^ is similar to that for E'^'^^. 



reasons of time-reversal symmetry it has to increase correspondingly before 
the jump to the next site. 

To analyse a possible time-dependence we recorded the time-dependence 
of the energy during residences by the ions. Here we restrict ourselves to 
the system LS— 2. A possible time-dependence would naturally depend on 
the typical residence time of an ion in a site. Therefore, we have grouped 
the different residence times into bins, namely S: r G [10 ps, 100 ps], M: 
r e [100 ps, 1 ns], and L: r G [1 ns, 10 ns]. The results are shown in Fig. |H 
For groups S and M we performed additional simulations with a sampling 
interval of 2 fs (during 1 ns). For group L the data from our main simulations 
is shown, with a sampling interval of 100 fs. Plotted times extend up to 10 % 
of the maximum duration, such that values are available from all residences 
in the group. 

The time dependencies of the energies from the start and from the end 
of the residence were both investigated. The positive direction of time was 
defined in both cases away from the nearest jump, such that the plotted time 
increases always towards the middle of the residence. Only the development 
from the start of the residences is shown here. In all cases the data from the 
end of the residences is identical within the errors, given by either error bars 
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Figure 4: Development of the energy of an ion directly after the jump into 
a site, at t = tjmp. Three sets of residences with different duration are 
distinguished, designated Short, Medium and Long. Behavior at the end of 
a residence before the jump is equivalent. 
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or symbol size. 

It is clearly seen that, if at all, E"^^^ and E^'^^ slightly change with time 
for the shortest residences. After their addition to the total energy E^°^ only 
unsystematic fluctuations remain. In conclusion no energetic relaxations are 
observed when an ion reaches a site, or before it leaves. The average site- 
energies are therefore the appropriate quantities for a deeper analysis of the 
behavior of the systems. 

4.3 Site— energy vs. thermodynamics 

The residence times are linked to the site-energies in an indirect way, via 
the height of the energy barrier between two sites. While it is probably 
not crucial, this complication is absent for the thermodynamic occupation 
probabilities of the sites, which should depend directly on their energies via 
a Fermi distribution. 

In Fig. [5] this prediction is tested for the system LS— 2. A Fermi function 
with fi = —8.4 eV is given by the solid line. The value of fi was chosen 
such that the predicted total occupation of all sites from the MD simula- 
tion corresponds to the number of alkali ions. However the occupation data 
for the individual sites from the MD simulation, given by the solid circles, 
deviate fundamentally from a Fermi function. No clear dependence of the 
relative occupation on the site-energy can be observed. For comparison, a 
system with the same number of sites and ions and with the site-energies 
E"*"* found in the MD system was propagated in a Monte Carlo simulation. 
The occupation data from this simulation, given by the open circles, nicely 
agrees with the theoretical Fermi function. 

Some minor deviations in the MD data might occur from the fact that in 
disagreement with the derivation of the Fermi relationship^ many sites may 
very rarely host two ions at the same time. However, as seen from Fig. [5]only 
4 sites have an average number of hosted ions which is significantly larger 
than unity. Therefore this effect cannot be responsible for the significant 
deviations. 

4.4 Site— energy vs. dynamics 

Another interesting observation in FigH] is the correlation of E'^^^ and E^^^ 
with the residence time. For long residences, represented by L, E'^^^ is smaller 
whereas E^'^^ displays the opposite trend. The data in Fig. H] also suggest 
that a minor correlation along this line seems to remain for Sites with 
short residence times have somewhat higher total energies E^"^. This has 
been systematically analyzed in Fig. [6] for the system LS— 2. For each site 
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Figure 5: Mean relative occupation of the sites in LS— 2 plotted against 
The solid line is a Fermi-Dirac function with fi = 8.4 eV, the open symbols 
were generated from a Monte Carlo model that recreates the MD system (see 
text). 
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the mean residence time Tresfor ions in a site is plotted against the network 
energy E^'^^ and against the site-energy E^°^. 

On average there is some increase of the residence time with increasing 
network energy. As directly indicated by the major scatter the correlation 
coefficient between logXres and i?'^^* is relatively small (0.28). The depen- 
dence of logTrea on E^°^ IS also present, but shows an even weaker correlation 
coefficient. (-0.18) A possible reason for these weak correlations is discussed 
below. 

5 Model energy landscape 

The deviations of ion behavior from expectations based on the average ener- 
gies, like the Fermi distribution for relative occupations, can be understood 
when considering the effect of ion-ion interactions. To visualize this effect we 
use a very simple model with sites and A^ — 1 ions. This choice reflects 
the observed small number of free sites. Let E'^'^^{s) denote the energy of an 
ion at site s due to the interaction with the network. Furthermore, i?™'^'*(s) 
denotes the additional energy of an ion at site s due to the interaction with 
the other ions under the condition that site t (t ^ s) is empty. One expects a 
small but significant dependence of _E'™'^'*(s) on t, e.g., via the distance of site 
t to site s. Furthermore, this term also contains possible correction effects of 
the network energy because the network structure may also depend on the 
actual ionic configuration. The energy ^vactotaij^^^, ^ Es^J^°''*(s) + ^™''*(s)] 
thus denotes the total energy of the system under the condition that site t 
is vacant. Correspondingly, the equilibrium probability p™'^'* that site t is 
empty is proportional to exp(— /5i?™*^'*°*'^^(t)). 

For this model one can ffist calculate the average energy E^°^{s) related 
to site s. It is given by 

E-(.) = E-%s) + ■ (1) 

Via the Boltzmann average this expression implies that the typical residence 
time is longer than the typical energy fluctuations via jumps of surround- 
ing particles. This ansatz is compatible with the simulation results because 
the fluctuations of the energy are much faster than typical residence times 
(t/zuc ^ Tres- The probability that site s is occupied by an ion can be calcu- 
lated as 

p->^(s) = 1 - J9™^'^ (2) 

In Fig. [7]p'°°(s) is compared to the Fermi distribution, based on the average 
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Figure 7: Comparison of the probability p'°°(s) for the ion population of 
site s in a simple vacancy model (see text) with its Fermi estimate. In one 
case (circles) the ion-ion interaction does not depend on the position of the 
vacancy, in the other case (squares) a weak dependence has been taken into 
account (see text for details) . 



site-energy E^°^{s) for the case that E'^'^^{s) is drawn from a box distribution 
[0,4A;sT]. Without the additional contribution of £'™'^'*(s) the Fermi distri- 
bution is indeed very well reproduced. Actually, for just a single vacancy the 
Fermi distribution is only an approximation but the difference to the true 
distribution is small. The situation dramatically changes when the addi- 
tional contributions E^^^'^{s) are taken into account. They are considered as 
random numbers, drawn from a narrower box distribution [0,0.5A;sT]. One 
can clearly see that now significant deviations from the Fermi distribution 
are observed. Thus, even a small interaction among the ions significantly 
invalidates the applicability of the Fermi distribution and renders averaged 
total site-energies unsuitable for the prediction of the population of that site. 

6 Discussion 

The site-energies that are the basis for this investigation of the energy land- 
scape for ionic transport can be determined with good precision as averages 
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from long molecular dynamics simulations. While there are fast fluctuations, 
the mean value of the site-energy is stable on the timescale of the typical 
residence of alkali ions in a site. The observed stability of the site-energies 
and the absence of relevant relaxations on longer time scales, i.e. on the scale 
of ionic diffusion or beyond, are necessary conditions for our analysis. 

A description of the energy landscape through a distribution of static 
site-energies, as it is used in most hopping models, is well justified by these 
results. Our analysis demonstrates that fundamental properties of the en- 
ergy landscape assumed by these models depend on the Coulomb interaction 
among the cations. Strong correlations in the network energies of neighbor 
sites are canceled out by the cation interaction to yield a nearly random 
landscape. Both the spatial correlations in E^'^^ and its anticorrelation with 
j^cat g^j^g reproduced by the counter ion model— by using full Coulomb in- 
teractions between mobile cations and fixed anions on a lattice. Arguably, 
any realistic energy landscape for a hopping model must at least implicitly 
incorporate effects of the cation interaction. One further example is the intro- 
duction of spatial correlations into the energy landscape of a model designed 
to reproduce the internal friction of mixed-alkali glasses.— 

Both observed correlations in the landscape can be tentatively explained 
by an argument similar to that used by Greaves in favor of conduction chan- 
nels.- He stated that any non-bridging oxygen (NBO) in the structure has to 
be shared between several cations in order to fulfill all coordination require- 
ments. Aided by the long range of the Coulomb interaction, the presence 
of an NBO can thus influence the energies of several neighboring sites alike. 
Higher numbers of available NBOs will similarly favor not a single site, but 
a region of the system. But such a region will be energetically favorable to 
the cations only just until their mutual repulsion balances the effect of the 
structure, giving rise to the observed anticorrelation. 

The correlations between neighboring sites and between the different con- 
tributions at one site are in turn essential for features of the landscape that 
directly affect the dynamics. While barrier heights for jumps between sites 
were not determined, the difference between two site-energies gives a lower 
bound for the barrier for the transition between them. Although the corre- 
lations in E^*^^ limit the occurring differences, the additional smoothing of 
the landscape by the anticorrelated E'^^^ is necessary to reach values as small 
as 0.33 eV for LS— 2. Interestingly, this value is significantly smaller than 
the macroscopic activation energy of 0.66 eV. Two reasons may play a role. 
One may expect that the barriers between the sites also make a significant 
contribution to the macroscopic activation energy. Moreover we have seen 
that the local dynamics is only weakly correlated with the local energy. 

In Fig. [6] the correlations between site-energies and residence times are 
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obscured by scatter that dominates the results for most sites. The larger 
values of observed for slow sites suggest a low density of non-bridging 
oxygen (NBOs) atoms and, correspondingly, a high density of bridging oxy- 
gen (BOs) nearby. This is fully compatible with the previous observation that 
sites with long residence times are surrounded by a larger number of BOs^. 
The physical reason is that the jump process of an ion is supplemented by an 
instantaneous local door-opening effect of oxygen atoms^i^. In case where 
the neighborhood contains BOs the door-opening effect is reduced, giving rise 
to a longer residence time. This structural trend is however only sufficient to 
indicate very slow sites, just as the network energies support only a general 
relationship, but no clear prediction of the dynamics at any individual site. 
The even weaker dependence of the residence times on E^°^ indicates that the 
average site energies do not capture all effects of the added cation interaction 
on the local dynamics. 

But also considering equilibrium properties the full site-energies, with 
the cation interaction included, have still proven insufficient to determine 
the behavior of ions at the individual sites. The most fundamental test 
for the importance of the site-energies in our systems is shown in Fig. 0, 
where the influence of the site-energies on their occupation probability is 
investigated. The relative occupation of the sites deviates fundamentally 
from the expected Fermi dependence. These observations clearly indicate 
that average single-site-energies, as defined here (thus, taking into account 
the long-range interaction) and used in many of the models of ion conduction, 
do not contain the complete information. 

The simple model analyzed in Fig. [7] suggests the reason: Coupling of the 
site-energies to other, distant parts of the system makes states of identical 
local occupation energetically non-equivalent. In consequence the statistics 
yielding e.g. the Fermi distribution are disrupted. In the glass simulated by 
the MD model, the fluctuations of the site-energy are fast compared to the 
duration of a residence. But for the case of only a few unoccupied sites the 
relaxation rate also strongly depends on the availability of free nearby sites. 
As jumps happen with delays much shorter than the total residence times 
once a free site is available, the effective timescale on which the site-energy 
acts on the occupation becomes comparable to the timescale of fast energy 
fluctuations. As a consequence only a weak correlation between residence 
time and site-energy is to be expected. 
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7 Conclusions 



We have shown that it is possible to describe the energy landscape for ion 
conduction in a glass in terms of stable average site-energies. The analysis 
of site-energies has established a central influence of the cation interaction 
on the basic properties of the energy landscape. Yet average site-energies 
including cation interactions are still insufficient to predict the ion dynamics 
directly. The dependence of the site-energies on the ionic configuration, 
and the low number of free sites, which regulates the possibility for jumps 
independently from the energies, are demonstrated as the likely causes. 

A direct description of the dynamics starting from site-energies will there- 
fore not be able to capture the correlations among the ions. Instead the falsi- 
fication of the single-energy picture suggests that alternative ways should be 
explored to build a description of the ion dynamics from information about 
the energy landscape. 

In particular a treatment in terms of vacancies seems promising— i^i^ 
because due to their small concentration interaction effects can be expected 
to be negligible. They offer an equivalent description of the ion dynamics, 
without loss of microscopic detail, because every jump of an ion can be 
replaced by a corresponding jump of a vacancy. Work along this line will be 
published elsewhere. 
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